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Abstract. We present and discuss isospectral quantum graphs which are not 
isometric. These graphs are the analogues of the isospectral domains in which were 
introduced recently in [1-5] all based on Sunada's construction of isospectral domains 
[H]. After presenting some of the properties of these graphs, we discuss a few examples 
which support the conjecture that by counting the nodal domains of the corresponding 
eigenfunctions one can resolve the isospectral ambiguity. 



1. Introduction 

M. Kac's classical paper "Can one hear the shape of a drum" triggered intensive 
research in two complementary aspects of this problem. On the one hand, a search for 
systems for which Kac's question is answered in the affirmative was conducted, and, 
on the other hand, various examples of pairs of systems which are isospectral but not 
isometric were identified. In the present paper we shall focus our attention to quantum 
graphs and in the following lines will review the subject of isospectrality in this limited 
context. The interested reader is referred to [1-11] for a broader view of the field where 
spectral inversion and its uniqueness are discussed. 

Spectral problems related to graphs emerge in two distinct ways. In the first, the 
spectrum of the connectivity (adjacency) matrix is considered. It represents a discrete 
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version of the Laplacian, and for finite graphs, the spectrum is finite. This set of 
problems is often referred to as combinatorial graphs. Quantum (metric) graphs are 
obtained by associating the standard metric to the bonds which connect the vertices. 
The Schrodinger operator consists of the one-dimensional Laplacians on the bonds 
complemented by appropriate boundary conditions at the vertices (see next section). 
The spectrum of quantum graphs is unbounded, and it displays many interesting features 
which made it a convenient paradigm in the study of quantum chaos |19j . 

Shortly after the appearance of Kac's paper, M. E. Fisher published his work "On 
hearing the shape of a drum" [T2j, where he addresses isospectrality for the discrete 
version of the Laplacian. Since then, the study of isospectral combinatorial graphs 
made very impressive progress. In particular, several methods to construct isospectral 
yet different graphs were proposed. A review of this problem can be found in |14j . 
In particular, a method which was originally put forward by Sunada [0] to construct 
isospectral Laplace-Beltrami operators on Riemann manifolds, was adapted for the 
corresponding problem in the context of combinatorial graphs. Here we shall go one 
step further, and show that it can be also adapted for quantum graphs. 

The conditions under which the spectral inversion of quantum graphs is unique 
were studied previously. In [151 [IB] it was shown that in general, the spectrum does 
not determine uniquely the length of the bonds and their connectivity. However, it 
was shown in [13] that quantum graphs whose bond lengths are rationally independent 
"can be heard" - that is - their spectra determine uniquely their connectivity matrices 
and their bond lengths. This fact follows from the existence of an exact trace formula 
for quantum graphs [18], |19]. Thus, isospectral pairs of non congruent graphs, must 
have rationally dependent bond lengths. The Sunada method, which is based on 
constructing the isospectral domains by concatenating several copies of a given building 
block, automatically provides us with graphs with rationally dependent lengths. An 
example of a pair of metrically distinct graphs which share the same spectrum was 
already discussed in [13]. In a previous report we have shown that all the known 
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isospectral domains in [SI H] have corresponding isospectral pairs of quantum graphs 
[2T] . Here, we shall take the subject one step further, and propose that isospectral 
graphs can be resolved by counting their nodal domains. That is, the nodal counts of 
eigenf unctions belonging to the same spectral value are not the same. The idea that 
nodal counts resolve isospectrality was suggested in [20] for a family of isospectral flat 
tori in 4-d, and was tested numerically. The present work offers both rigorous and 
numerical evidence to substantiate the validity of this conjecture in a few examples of 
isospectral graphs. For these examples the nodal counts differences ocur on a substantial 
fraction of the spectrum. 

The paper is organized in the following way. For the sake of completeness we shall 
give a short review of some elementary definitions and facts on quantum graphs. We 
shall then show how pairs of isospectral domains in can be reduced to isospectral 
pairs of quantum graphs, and discuss their spectra and eigenfunctions. Finally, we 
shall study the nodal domains of these eigenfunctions and show that by counting nodal 
domains one can resolve the isospectral ambiguity of the graphs presented. 

1.1. A short introduction to quantum graphs 

We consider finite graphs consisting of V vertices connected by B bonds. The V x V 
connectivity matrix will be denoted by Cij : Cij = r when the vertices i and j are 
connected by r bonds, and it vanishes otherwise. The group of bonds which emerge 
from the vertex i form a "star" which will be denoted by S^'^\ The valency Vi of a vertex 
is defined as the cardinality of the star S**-*-* and Vi = '^jCij. Vertices with Vi = 1 
belong to the graph boundary. The vertices with Vi > 1 belong to the graph interior. 
The bonds are endowed with the standard metric, and the coordinates along the bonds 
b are denoted by Xf,. The length of the bonds will be denoted by Lb, and the total length 
of the graph is C = J2b^b- 

The domain of the Schrodinger operator on the graph is the space of functions which 
belong to Sobolev space H^{b) on each bond b and at the vertices they are continuous 
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and obey boundary conditions as is mentioned in ([T]). The operator is constructed in 
the following way. On the bonds, it is identified as the Laplacian in 1-d — It is 
supplemented by boundary conditions on the vertices which ensure that the resulting 
operator is self adjoint. We shall consider in this paper the Neumann and Dirichlet 
boundary conditions: 



Neumann : - — ^bixb] 



= 

Xi,=0 



Dirichlet Vi : '^/'^(xfe)!^^^^ = . (1) 

The derivatives in ([T]) are directed out of the vertex i. Comment: The Neumann 
boundary conditions will be assumed throughout, unless otherwise stated. A wave 
function with a wave number k can be written as 

'^h{xh) = . \ J {4>i sin k{Lb - Xb) + 4>j sin kxb) (2) 
sm fhi^b 

where b connects the vertices i and j, where the wave function ipf, takes the values 
(pi and (f)j respectively. The form ([21) ensures continuity. The spectrum {kn} and the 
corresponding eigenfunctions are determined by substituting ([2]) in ([T]). The resulting 
homogeneous linear equations for the (pi are written as 

B 

Wl<i<V : = , (3) 

and a non trivial solution exists when 

/(Li, ■■■,LB;k) = det A{Li, ■ ■ ■ , Lb; k) = . (4) 

The spectrum {k^}, which is a discrete, positive and unbounded sequence is the zero 
set of the secular function f{Li, ■ ■ ■ ,Lb', k). The secular functions of the type (jl]) have 
poles on the real k axis, which renders them rather inconvenient for numerical studies. 
The secular function can be easily regularized in various ways, (See e.g., [THl 1^). 

It is easy to show that the complete wave function can be written down in terms 
of the vertex wave functions at the interior vertices with Vi > 3 only. In the sequel we 
shall denote their number by Vint- This reduces the dimension of the matrix A above 
from V to Vint- 
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The nodal domains of the eigenf unctions (the connected domains where the wave 
function is of constant sign), are of two types. The ones that are confined to a single 
bond are rather trivial. Their length is exactly half a wavelength and their number is 
on average The nodal domains which extend over several bonds emanating from a 
single vertex vary in length and their existence is the reason why counting nodal domains 
on graphs is not a trivial task. The number of nodal domains in a general graph can be 
written as 

= ^ E E { L^J + ^ (l - (-l)L'^Jsign[0.]sign[0,]) \-B + V. (5) 

where [xj stands for the largest integer which is smaller than x, and (pi, (pj are the values 
of the eigenfunction at the vertices connected by the bond b [22J. holds for the case 
of an eigenfunction which does not vanish on any vertex: (pi 7^ 0. 

Recently Schapotschnikow [23j proved that Sturm's Oscillation Theorem extends 
to finite tree (loop-less) graphs: the number of nodal domains of the n'th eigenfunction 
(ordered by increasing eigenvalues) is n. Berkolaiko |;24j have shown that the number 
of nodal domains is bounded to the interval [n — l,n] where I is the minimal number of 
bonds which should be cut so that the resulting graph is a tree. 

Nodal domains can be also defined and counted in an alternative way which makes 
use of the vertex wave functions {(pi} (see ([2])) exclusively: 

A nodal domain consists of a maximal set of connected interior vertices (vi > 3) 
where the vertex wave functions have the same sign. 

This definition has to be modified if any of the (pi vanishes. Then, the sign attributed 
to it is chosen to maximize the number of nodal domains [25]. 

We thus have two independent ways to define and count nodal domains. To 
distinguish between them we shall refer to the first as metric nodal domains, and the 
number of metric domains in the n'th eigenfunction will be denoted by z/„. Berkolaiko's 
theorem states that (ri — /) < z/„ < n. The domains defined in terms of the vertex wave 
functions will be referred to as the discrete nodal domains. The number of discrete 
nodal domains of the n'th vertex wave function will be denoted by The sequences of 
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metric and discrete nodal domains counts {z/„} and {/in} are the main objects of study 
of the present paper. 

2. Isospectral quantum graphs 

The first pair of isospectral planar domains which was introduced by Gordon, Web and 
Wolpert pLj is a member of a much larger set which was discussed in [3l H] . This was 
extended in [5] to include domains which differ in the distribution of boundary conditions 
(Dirichlet or Neumann) along their boundaries. The common feature of these sets of 
pairs of isospectral domains is that they are constructed using the Sunada method [6], 
and they share a few important and distinctive attributes: 

• The domains are constructed by concatenating an elementary "building block" in 
two different prescribed ways to form the two domains. A building block is joined to 
another by reflecting along the common boundary. The shape of the building block 
is constrained only by symmetry requirements, but otherwise it is quite general. 

• The eigenfunctions corresponding to the same eigenvalue are related to each other 
by a transplantation. That is, the eigenfunction in a building block of one domain 
can be expressed as a linear combination of the eigenfunction in several building 
blocks in the other domain. The transplantation matrix is independent of the 
considered eigenvalue. 

• The construction of these pairs reflects an abstract algebraic structure which was 
identified by Sunada [6]. 

An example of an isospectral pair of domains in and its building block is given 
in figure ([11(a)). This is a pair whose construction is described in [3]. It is denoted there 
by 73. Other examples are displayed in e.g., [21 El [5]. 

We can construct metric graphs which are analogous to these domains, by replacing 
the building blocks by appropriate graphs which preserve the requires symmetry. As an 
example, the triangular building block in figure (Pl^a)) can be replaced by a 3 - star with 
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Figure 1. (a) Planar isospectral domains of the 73 type, (b) Reducing the building 
block to a 3-star.(c) The resulting isospectral quantum graphs. 

bonds of lengths a, b and c as shown in figure (mb)). This yields the pair of isospectral 
but non isometric graphs shown in figure (mc)). (In drawing figure dUc)) we took 
advantage of the fact that the "angles" between the bonds are insignificant). Note that 
the two graphs share the same connectivity matrix (they are topologically congruent) 
however they are not isometric. As a matter of fact the right and the left graphs are 
interchanged when the bonds "6" with "c" are switched. Thus, the lengths b and c 
must be different to ensure that the two graphs are not isometric. This is an example 
of an asymmetry requirement which the building block must satisfy. All other pairs of 
isospectral domains proposed in [3j can be converted to their analogous quantum graphs 
(see [2l] for details). Two examples are shown in figure (Ej). 

The most simple pair of isospectral domains can be built by applying the Sunada 
method to the Dihedral group D4. The resulting graphs are shown in figure ([3]), where 
the letters D or N specify the boundary conditions at the boundary vertices. (The 
boundary conditions at the interior vertices are always Neumann). Note that the 
two graphs are connected differently, in contrast with the previous example where the 



topology of the two graphs is the same. In [Appendix A| we give a brief description of 
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Figure 2. The isospectral pairs 183 (left) and ISs (right) [3] 
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the construction of this pair, as an example of the Sunada method. 

In a way of illustration we shall write explicitly the secular functions for the graphs 
shown above. Starting with the 7s pair, we have Vint = 7 (the interior vertices are 
marked by a black dot and are enumerated in figure 1(c)). The secular function for the 
graph / is written in terms of the matrix 
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Here, 



a{a; k) 



sm{2ka) 



Pib; k) 



sm.{2kh) 



7(a; k) 



sin(2A;c) 



(7) 



and 



^(a, 6, c; /c) = tan(aA;) + taxiipk) + taxi{ck) . (8) 

The corresponding vertex eigenfunction cj) ^ = {4>i , ■ ■ ■ , (p-j^) is the eigenvector of 
A{a, b, c; with a vanishing eigenvalue. The A matrix for the graph // is obtained 
from the matrix A{a, c, b; kn) by interchanging c and b. Exphcit computation shows that 
det A{a, b, c; k) = det A{a, c, b; k). This proves the isospectrahty of the two graphs. 

The sum of the elements of any of the columns of A{a, b, c; k) is ^. Hence the vector 
1 = (1, 1, 1, 1, 1, 1, 1) is an eigenvector with an eigenvalue ^. The values of k for which 
^{a, b, c;k) = are in the spectrum of both graphs. They correspond to eigenf unctions 
which are the same on each of the 3-stars - the building block of the complete graphs 
(figure 1(b)). As a matter of fact, the condition ^(a, b, c;k) = is the secular equation 
for this 3-star with bond lengths a,b,c [I9]. This subset exhausts 1/7 of the spectrum 
of the graphs, and in this set the transplantation property is trivial. 

The transplantation property which is basic to the proof of isospectrality for the 
domains in [3] can be explicitly formulated by the transplantation matrix 



/ 1 1 1 

1 1 1 

10 10 1 

1 1 1 
10 10 10 
1 1 1 

y 1 1 1 



\ 



so that (p" = T(j)^ . 



(9) 



We observe that the inverse transplantation is also effected by T: 



II 



even 



though T is not self inversive. The fact that T induces the transplantation in the 
two directions, implies that the vertex wave functions (f) must be eigenvectors of 
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(as long as the eigenvalue is not degenerate). The wave functions are defined up to 
normalization, and therefore the corresponding eigenvalues can be different than unity. 
The spectrum of consists of the eigenvalue 9 and the six-fold degenerate eigenvalue 
2. Thus, the vertex wave functions are either proportional to 1, or belong to the 6 
dimensional subspace of vectors orthogonal to 1. This observation will be used in the 
next section when the nodal structure of the eigenf unctions is to be discussed. 

Other examples of pairs of isospectral graphs such as e.g., the pair described in [13] 
can be obtained from the isospectral graphs constructed above, by setting the length of 
the c bond to zero. This was explicitly shown in [21]. Using similar special cases one 
can generate a rich variety of isospectral but not isometric quantum graphs. 

Finally, it should be emphasized, that the skeleton of interior vertices form graphs 
which are identical to the topological, colored graphs which were introduced in [1] to 
express the transplantation properties of the isospectral domains in M^. The quantum, 
metric graphs can be obtained simply by completing each topological vertex to a 3- 
vertex star graph with lengths a, b, c. Topological bonds are associated with lengths 
according to their colors, and bonds are added to the topological vertices with valency 
less than three, to complete them to 3- stars. This proves the statement made above 
that the computation carried out here can be repeated for any of the topological graphs 
shown in [1]. The graphs corresponding to the isospectral domains 133 and 138 are 
shown in figure [21 




N 



I II 



Figure 3. The isospectral pair with boundary conditions. D stands for Dirichlct and 
N for Neumann. 
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The Dihedral graphs are much simpler, with Vint = 2. Their secular determinants 
can be shown to be identical functions of k. This is an explicit demonstration that the 



graphs are indeed isospectral. The secular equation is obtained from 
f n — (3 + "y) —a ^ 

A{a, b, c; k) = 



-a ?7 + (/5 + 7) 
(3ib;k) 



"^""'^^ sin(2A;a) ' "^"''^^ sin(2A;6) 
7]{a, b, c; k) = cot{2ka) + cot(2A;6) + cot(2A;c 



7(a; k) 



sin(2/cc) 



(10) 



As can be easily seen, det A has poles, and to get a regular secular function we have to 
multiply it by sin(2fca) sin(2A;6) sin(2A;c). It takes the form 



/(a, b, c; k) = sm{2ka) (-2 + 2 cos(2A;6) cos(2A;c) - 3 sin(2A;6) sin(2A;c)) 

+ 2 cos(2fca) sin(2A;6 + 2kc) 



:iii 



The transplantation matrix for this pair of graphs is derived in Appendix A It 



reads 



T 



V2 



\ 



so that 



T0 ^ . 



1 -1 

V ^ W 

The eigenvectors ^'^^ point at the directions 6 with tan(^ ^'^^) 



(12) 



The 



transplantation ( |T2i) implies that is obtained from ^ by a rotation of | 
counterclockwise. Direct substitutions shows that 



tan(6' \kn)) = g{a,b,c; /c„) 



(13) 



/ , ,x X . , /cos(2A;6) — 1 cos(2A;c) — T, , , 
g{a, b, c; k) = cos(2fca) + sm(2fca) . ; \ . / . (14) 



where, 



sin(2/c6) sin(2/cc) 
The explicit form of the functions f{a,b,c;k) and g{a,b,c;k) will be used in the 
discussion of nodal counting in the next section. 
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3. Nodal counts and the resolution of isospectrality 

The real eigenfunctions of Laplace-Beltrami operators on two dimensional manifolds 
display usually an intricate sets of nodal domains - which are the connected domains 
where the sign of the eigenfunction is constant [27]. Sturm's oscillation theorem for 
systems in 1-d, and its extension by Courant to any dimension, establish the connection 
between the number of nodal domains and the spectrum: The number of nodal domains 
Un of the n'th eigenfunction is bounded by n. (the eigenfunctions are arranged by 
increasing value of their eigenvalues). Courant's theorem for combinatorial graphs and 
for quantum graphs were proved in [22] , respectively. 

Given a pair of isospectral domains. Are the sequences of nodal counts {z/„}^^ 
identical? In other words, can one use the information stored in the nodal sequences to 
resolve isospectrality? This question was recently discussed in the context of isospectral 
flat tori in R*^ with d > 4 [20], and numerical as well as analytical evidence was brought 
to substantiate the conjecture that the nodal sequences resolve isospectrality. In this 
section we shall provide rigorous as well as numerical evidence to show that the same 
is true for isospectral quantum graphs of the types discussed above. We shall start by 
discussing the discrete nodal sequences and then proceed to the metric nodal sequences. 

3.1. The discrete nodal sequences 

The discrete nodal domains were deflned as the maximally connected sets of interior 
vertices with vertex wave functions of equal sign. Their number is denoted by 
and the nodal sequence is We shall prove in the next subsection that for the 

dihedral graphs half of the entries in the sequences of discrete nodal counts are different. 
We are not able to provide a similar proof for the more complex graphs, but we shall 
bring numerical evidence which supports the conjecture that their sequences of discrete 
nodal counts are different. 
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3.1.1. The dihedral graphs The number of discrete nodal domains for this pair of 
dihedral graphs (see figure [3]) can be either one or two, depending on whether the signs 
of the components of the (two dimensional) vertex eigenfunctions have the same sign or 
not. In other words (see ( |T3|) ). the number of nodal domains depend on the quadrant 
in the (0i,02) plane where the eigenvectors point: /i = 1 in the first and the third 
quadrants, and /i = 2 in the second and the fourth quadrants: 

/i„ = 1 + ^ (1 - sign(tan6'„)) . (15) 

The transplantation implies that the eigenvectors 0/^ are obtained by rotating 0/ by | 
counterclockwise. Therefore, /x^ ^ fi^J if the transplantation rotates the vectors across 
the quadrant borders. In other words, 

f^i + l^n ^ tan(^0 e {(-1, 0) U (1, oo)} (16) 

This observation is essential to our discussion since it expresses the problem of nodal 
counting in geometrical terms. 

It is convenient to construct finite subsequences {fi^} and {fJ^n} of discrete nodal 
count of graphs I and II, restricted to the spectral points in the interval < kn < K. 
We denote the number of terms by N{K). Define 

^^^^ = NlK) ^^^^ ■ ^ -""'^ • ^^^^ 

We shall now prove 

Theorem 1 Consider the dihedral graphs I, II discussed above with rationally 
independent bond lengths a, b, c. Then 

hm PiK) = i. (18) 

K^oo Z 

Proof of theorem [I].- The rational independence of a, b, c implies that the eigenfunctions 
never vanish on the inner vertices and hence that the spectrum is simple. 

In order to study P{K) above, we consider the distribution of the directions of the 
eigenfunctions 0^ in the spectral interval. Using flllfl3[ [H]) we get 

S{x — g{a,b,c; k)) (19) 



h{x; K) = {6{x - tanen))K = -tttttt / d/c 6{f{a, b, c; k)) 

J^{1^) Jo 



df 



dk 
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In taking the limit ^ oo we use N(K) = 2a+^+2c^ jjgj^ Moreover, since a,b,c 
are assumed to be rationally independent, k creates an ergodic flow on the 3-torus T3 
spanned by r = 2ka mod27r, s = 2kb mod27r, t = 2kc mod27r Ergodicity implies 
that the integral over k in ( fT9|) may be replaced by an integral over T3 leading to 



h{x) 



2{a + b + c)n~ jo 



1 



(*27T p27V nliT 

dr / ds / dt S{f{r,s,t)) 
Jo Jo 

ds dti drS{fir,s,t)) 
2{a + b + c)ti^ Jo Jo \Jo 

+ I drSif{r,s,t)) 



df 



dk 
df 



dk 



S{x - g{r,s,t)) 
6{x - g{r,s,t)) 



df 



dk 



6{x - g{r,s,t)) 



i-2n r'^n 

2{a + b 

Now we note that under the transformation r 



2 i'2n i'2'K 

— — — / ds dt[Ii{s,t;x) + l2{s,t;x)] 
b + cjvr^ Jo Jo 



(20) 



(r + 7t) mod27r we have 



fir',s,t) = -f{r,s,t) 

df , , , d/ , , 
dfc^^'^'') = -d^^^'^'') 

g{r',s,t) = -g{r,s,t) 



Thus, we conclude that 



h(s,t; x) = his.t] -x) 



h{x) = h{—x) 



-1 



h{x)dx + / h{x)dx 



h{x)dx + / h{x)dx = - 
-1 Ji 2 



lim P{K) = - 

K^co 2 



□ 



3.1.2. The graphs Ts, 183 and ISs At this stage we are not able to prove the validity 
of the conjecture that the nodal sequences resolve the isospectrality of the graphs Ta, 
183 and ISg. However, numerical tests show that for an appreciable fraction of the wave 
functions, the nodal counts are distinct. We shall start the discussion by considering 
the case 73, for which we can turn the counting problem into a geometrical problem. 

The vertex wave functions which correspond to the same eigenvalues are related by 
transplantation: 

$'' = T$' ■ 0^ = T0^^ (21) 
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— * 

where (p stand now for the vertex wave function restricted to the 7 interior vertices in 
these graphs. We have already shown that eigenvectors (f)^'^^ are either proportional to 

1 (the vector with constant entries)- in which case both have a single nodal domain, 
and = 0, or the are orthogonal to 1. The 6 eigenvectors of with eigenvalue 

2 span the orthogonal subspace: three of them correspond to the eigenvalue +\/2 of T, 
and the other three correspond to the eigenvalue — -\/2 of T. Any vertex wave function 
which is not proportional to 1, can be written as 

^ = cosa 1+) + sina |-) (22) 

where |±) stand for arbitrary normalized vectors in the three dimensional subspaces 
mention above, and a e [0, 27r] . Thus, 

= r0 = y2(cosa 1+) -sina |-)) (23) 

and the difference in nodal numbers comes only through the change of sign of the linear 
combination. We are not able to make further progress beyond this point, and we shall 
therefore summarize the numerical findings. 

The graph 7^ - Based on the study of approximately 6600 eigenfunctions we find 
that the difference between the discrete nodal counts can take the values 0, ±2 only. 
The nodal counts are different for ~ 19 % of the wave-functions. 

The graph 183 - Based on the study of approximately 2000 eigenfunctions we find 
that the difference between the discrete nodal counts can take the values 0, ±2 only. 
The nodal counts are different for ~ 22% of the wave-functions. 

The graph ISg - Based on the study of approximately 4700 eigenfunctions we find 
that the difference between the discrete nodal counts can take the values 0, ±2, ±4 only. 
The nodal counts are different for ~ 22 % of the wave-functions. 

This result is certainly encouraging but not sufficient, and we are trying various 
options to substantiate it more rigourously. 
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3.2. The metric nodal sequences 

The discussion of the metric nodal sequences is facihtated in the present examples 
of isospectral graphs since most of them are trees. Then, by making use of 
Schapotschnikow's theorem [23], the metric nodal counts Un = n. This applies in 
particular to the isospectral graphs of the types 7 3, 183 and ISg which are tree graphs and 
~ ~ ^- The metric counts do not resolve the isospectrality in these case. However, 
less trivially connected graphs can be constructed by replacing the basic building block. 
For such graphs Schapotschnikow's theorem does not apply, and we expect their metric 
sequences to be different. Below we shall show that this is true for the pair of dihedral 
graph where graph II is not a tree. 

Consider the pair of dihedral graphs, and let {i^^^} denote the metric nodal counts 
of graphs I and and II. Using ([5]) we express the difference — u^'^^ as 

= vi- vl! (24) 
= ^ [ 1 - sign(02,nSin(2A;„a)) + sign(02ij sin(2A;„6)) + sign(02f„ sin(2A;„c)) ] 

To obtain ([21]) a slight modification of ^ was needed: the term in the sum that 
corresponds to a bond with boundary vertex that has Dirichlet boundary conditions 
should be modified to be instead of + | - (-l)^^^sign[0j]sign[0 

In addition while deriving f[2^ we made use of the identity [^J — 2[|J = ^(1 — 
sign(sin(2a;))) and of the freedom to choose the sign of the first component of the wave- 
functions to be positive. A theorem by Berkolaiko [2^ guarantees that (5z/„ can only 
take the values and 1. We would like to compute the distribution of these values on 
the spectrum. For this purpose we consider the spectral interval < /c„ < K and study 
the function 

Q{^\ K) = ]^tl{n < Nik) : Su^ = x} (25) 
and prove the following theorem. 

Theorem 2 Consider the dihedral graphs I, II discussed above with rationally 
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independent bond lengths a,b,c. Then 

3 lim Qix,K) = Q{x) and g(0) = Q(l) = i. (26) 

K^oo Z 

Proof of theorem\^ In order to study the Q{x,K) above we consider the distribution 
of (5z/„, which is given by the integral 

df 



h{x;K) = {6{x-5u^)) 



K 



K 



dk 



6{x- 6iy{a,b,c;k)) , (27) 



dk6{f{a, b, c; k)) 

where the spectral secular function /(a, b, c; k) is defined in (iTTl) above, and 5i/(a, b, c; k) 
coincides with 5z/„ for k = k^ and can be written explicitly as 

Ma,&,c;A;) = i[l (28) 



— sign ( cos(2A;a) sin(2/ca) + sin^(2A;a)(cot(2/c6) + cot(2/cc) — 
sin^(2A;6) sin(2A;c) / 1 



+ sign 



+ sign 



sin(2A;6) + sin(2A;c) Ysin(2A;a) 
sin(2A;6) siii^{2kc) 



sin(2A;6) sin(2A;c)' 
+ cot(2A;a) + cot(2A;6) + cot(2A;c) 



+ cot(2fca) + cot(2fc6) + cot(2fcc) 



sin(2A;6) + sin(2/i;c) \sin(2A;a) 
Following the same route as in the proof of theorem ([T]) we take the limit K ^ oo while 

making use of the ergodic theorem, and replace the k integration by an integration over 

the 3 - torus with coordinates r = 2ka mod27r, s = 2kb mod2'7r, and t = 2kc mod2'7r. 



h{x) 



2a + 26 + 2c7r3 



2n 



dr 



2-K 



ds 



2n 



df 



dk 



6{x — 5z/(r, s, t)) 



dM(/(r,s,t)) 

= Qi-l)5ix - (-1)) + Qi0)5{x) + Q{l)5{x - 1) + Q(2)5(a; - 2) (29) 

The last line follows from the fact that x is an integer which is written as half 
the sum of four unimodular numbers. Next we note that under the transformation 
r I— >• r' = (— t) mod27r, s ^ s' = (— s) mod27r, t ^ t' = {—t) mod27r we have 

f{r',s',t') = -/(r,s,t) 

df . , , df . , 

-(r,.,t) = -(r,.,t) 

6h'{r' , s', t') = 1 — 5z/(r, s, t) 

Thus, we conclude that Q{—1) = Q{2) A <5(0) = (5(1). But due to Berkolaiko's 
theorem [21] Q(-l) = Q{2) = 0. Therefore, 



Q{0)=Q{1) = - 



□ 
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4. Discussion and summary 

There is now a growing amount of evidence that nodal-count sequences store information 
on the geometry of the system under study, which is similar but not equivalent to the 
information stored in the spectral sequence. In a recent paper [2H] it was shown that 
the nodal count sequence for separable Laplace - Beltrami operators can be expressed in 
terms of a trace formula which consists of a smooth (Weyl - like) part, and an oscillatory 
part. The smooth part depends on constants which can be derived from the geometry 
of the domain, and the oscillatory part depends on the classical periodic orbits, much 
in the same way as the spectral trace formula. However, what we have shown in the 
present paper is that the information stored in the two sequences is not identical: For 
the isospectral pairs of graphs considered here, the nodal-count sequences are different 
in a substantial way. In this respect, the nodal sequence resolves isospectrality. 

Till this work was done, the conjecture that isospectrality is resolved by counting 
nodal domains was substantiated by numerical studies only. Here we presented for the 
first time a system where this fact is proved rigorously. The main breakthrough which 
enabled the proof was by formulating the counting problem in a geometrical setting. We 
hope that this way will pave the way to further analytical studies, where more complex 
systems will be dealt with. 

Finally, we would like to mention a set of open problems which naturally arise in the 
present context: Can one find metrically different domains where the Laplacians have 
different spectra but the nodal counting sequences are the same? A positive answer is 
provided for domains in 1 dimension (Sturm) or for tree graphs (Schapotschnikow). Are 
there other less trivial examples? 
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Appendix A. Construction of the isospectral dihedral graphs 




Figure Al. A square with its 4 possible axis of reflection 

The Dihedral group D4 is the symmetry group of a square. Let a be the 
transformation of rotating the square n/2 counterclockwise and b be the transformation 
of reflecting it along the x axis. The Dihedral group is generated by a, b which obey 
= 1, 6^ = 1, bab = a~^. The 8 group elements of D4 are {1, a, a^, a^, = 
b, Vy = aba~^, r„ = ab, = ba}. r^, Vy, r„, are reflections along the axis 
X, y, u, V respectively (figure lAll) . Note that the multiplication should be read from 
right to left since the elements are transformations. This group has 5 irreducible 
representations (irreps). Four of them are one dimensional and one, which is of our 
interest, is two dimensional. The two dimensional representation is: 
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1 ^ 



1 
1 

-1 

1 

1 

1 , 



V 



\ 



1 
-1 

-1 

ly 

-1 
-1 



-1 







-1 



1 
-IJ 



J 



(A.l) 



Denote the appropriate basis by which the representation looks as above by 



= {^{, \l/2}- By change of basis using the matrix T = 



we obtain the similar representation 



V 



-A 



[X 



T-^xT) 



1 ^ 



( 



( 



1 
1 

-1 

1 

1 
-1 



( 



1 

-1 

1 

1 

-1 
1 







' -1 ^ 







(A.2) 



Denote the second basis by B = {\E'{ , \E'2 }• So that 



T 



Notice that in the first form of the representation and Vy are diagonal and in the 
second one r„ and r„ are diagonal. This will be exploited soon. 

The graph in figure IA2( a) obeys the dihedral symmetry. Examine the 
eigenfunctions of the Schrodinger operator with Neumann boundary conditions on all 
the vertices. The set of all the eigenfuncions of a certain eigenvalue A form a finite 
dimensional vector space V with some dimension dimV = n. This vector space is 
invariant under the action of the dihedral group and therefore V is the carrier space 
of some representation of D4. We are interested in the case that this representation is 
the two dimensional irrep of D4 or that this representation is a reducible representation 
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J H 



(b) 

Figure A2. A graph which obeys the dihedral symmetry. In (a) the lengths of the 
bonds are marked. In (b) The points on which we can deduce boundary conditions are 
marked. 

that contains the two dimensional irrep. Consider the basis = "^2} of this 2d 
irrep (note that now have the meaning of wave functions on the graph in figure 

IA2[) . From the matrix which corresponds to in this basis (see ( lA.ll) ) we deduce that 
vanish on the axis x (as an odd function with respect to that axis) and as for \l/2, its 
derivative vanishes on the x axis (since \l/2 is an even function with respect to that axis 
and therefore its derivative is odd). Using the notation from figure R2] (b). we have: 



E,F,K,L 



dx 



(A.3) 



E,F,K,L 



Similarly, examining the matrix of Ty in this basis flA.ip yields 

= =0 



B,C,I,H 



(A.4) 



We can get similar observations by looking on the matrices of r^, r^, this time for the 
basis B^^ = {^{, ^^^} ([Q) . We get the following 
d 



dx 



II 







D,J 



2 Id, J 







II \ 



\A,G 



dx 



II 



(A.5) 
(A.6) 



A,G 
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Now, a pair of isospectral graphs will be constructed out of the complete graph. 
Each of the graphs in the pair will be a subgraph of the complete graph. Graph I is the 
graph that consists of the vertices B,C,D,E,F and all the bonds which connect them. 
Graph II is the graph that consists of the vertices A,B,C,D and all the bonds which 
connect them (figure EH]) . 

We can define the subgraphs in another way which will become convenient later: 
using the axis x, y, u, v as in figure ED we define two sets of coordinates (x, y) and {u, v) 
and use them to examine the complete graph. Graph I is obtained when restricting 
ourselves to the region Dj = {x ^ A y ^ 0} of the complete graph. Graph 
II is obtained when restricting to the region Djj = {u ^ A v ^ 0}. For a 
certain eigenvalue A of the complete graph whose representation contains the 2d irrep 
discussed before, we can define appropriate eigenf unctions on the two subgraphs by 
ijj^ = ,ip^^ = ■ The boundary conditions of ip^^il^^^ are obtained from 

f[0|) . f[01) . f[0|) . f[0|) and are shown in figure |M1 




2b 



D a I I a N 




2c 

I II 

Figure A3. The isospectral pair with boundary conditions. D stands for Dirichlct 
and N for Neumann. 

So far we have established the ground for the proof of the following: 

Theorem 3 Let a' denote the spectra of graph I ( with its boundary conditions ) and a'^ 
denote the spectra of graph II. Then with account of multiplicities. 

Proof of theorem We will make use of the results above in order to construct a 
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transplantation which proves the theorem. Let A G cr^. Denote by the space of 
eigenfunctions which belong to the eigenvalue A. V?/'^ G define the functions \&f , \&2 
on the complete graph by the following: 

ip^x.y) for X ^ A y^O 

—ih^ix, —y) for X ^ A y ^ 

ip\—x,y) for X ^ A y^O 

—tp\—x, —y) for X ^ A 2/^0 
mx,y) = -^[{y,-x) (A.8) 

Where we have used again the coordinates (x, y) as defined trivially by the axis x, y 
(scaling is obviously not important for our purposes). It can be verified that \E'f, 
'^2 valid eigenfunctions of the eigenvalue A defined on the complete graph. If 
\E'((x, y) = c\I'2(x, y) for some c G M then we get that ip^i^x, y) = cip^y, x) ^ ip^ = (this 
can simply be concluded by considering the boundary conditions). Therefore \I'2 
are independent eigenfunctions and form a basis = {\E'(, \I'2} for the appropriate 
eigenspace. This is exactly the basis of the first form of the 2d irrep of D4 (see lA.ll) . 



(This can be verified for example by examining the representation of the generators a, b 



of D4). Now, we construct the basis B = {\E'{ , \E'2 } by 



1 

V2 



-1 



V 



V*2 y 



. And now define the function ip^^ on graph II by ip^^ = 



T 



We get that tp^^ is an eigenfunction of graph II which belongs to the eigenvalue A and 
obeys the appropriate boundary conditions ^ A G ct/j. 

Following this procedure, one can verify that T is the transplantation matrix 
of the vertex wave functions. We can similarly describe the inverse transplantation. 
The inverse transplantation is linear and therefore if we have an independent set of 
eigenfunctions (with the same eigenvalue) they would remain independent after the 
transplantation. Thus, the degeneracies of A in aj and ajj are equal. The existence of 
the inverse transplantation also proves that VA G an ^ A G a/ and the degeneracy is 
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again the same. Therefore ar = cru with account of multiphcities. □ 
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